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Abstract 

We investigate tlie two-loop gap equation for the thermal mass of hot massless 
g'^(p'^ theory and find that the gap equation itself has a non-zero finite imaginary 
part. This indicates that it is not possible to find the real thermal mass as a solution 
of the gap equation beyond order in perturbation theory. We have solved the 
gap equation and obtain the real and the imaginary part of the thermal mass which 
are correct up to g"^ order in perturbation theory. 
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It is well-known [||] that if a theory contains massless bosonic field such as, QCD 
or massless scalar theory with (7^0^ interaction, then at very high temperature (T) the 
perturbative computations beyond certain order of coupling constant are afflicted with 
infrared (IR) singularities. In the case of massless (7^0^ theory, the one loop contribution 
to two point function shows that the fields are screened and the screening mass (Debye 
mass) is found to be of order gT p|, ^. However, the result of two loop corrections to 
it is found to be IR divergent. A natural way of avoiding this IR divergences in the two 
loop computation is to use the dynamically generated one loop thermal mass gT as the 
lower cut-off of the momentum integration. As a result one finds the appearance of a 
new g^ order correction to two point function which, although consistent with the spirit 
of perturbation theory, has not been predicted from the usual perturbative expansion in 
powers of coupling constant {g) at zero temperature. In addition to that there are infinite 
number of higher order diagrams that contribute to this particular g^ order which in turn 
is the signature of the break-down of usual perturbation theory at very high temperature. 
Moreover, it also suggests that one has to resum this infinite number of diagrams to 
correctly calculate this g^ order contribution 0, ^, ^ to two-point function. 

However remaining within the framework of perturbation theory one can in principle 
be able to calculate the thermal mass to any order of coupling constant by solving the gap 
equation @, 0. In order to obtain the gap equation, the functional integral formulation 
may be used, and we shall briefiy discuss this method following Jackiw et. al. Consider 
a massless g'^cj)'^ theory in d dimension (d = 4 — 2e) described by the lagrangian 

c = U,<p^^<p-^^'^^^ct>\ (1) 

and the partition function is given by 

Z = j r'0e*^('^) , (2) 

where 5(0) = / d'^xC{(f){x)). We introduce a loop counting parameter / and write down 
Z as 

Z = J P0et^(^*) . (3) 

In the usual perturbation theory we separate the quadratic part of 5(0) and expand the 
exponential of the remainder in powers of Z. To obtain a gap equation for a possible mass 
m, we add and subtract Sm = f d^x(f)'^{x), which of course changes nothing (at least 
in the classical level, the equation of motion remains the same). 

S = S + S,n-Sm (4) 
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We recognize the loop expansion by expanding S + Sm in the usual way, but taking —3^ 
as contributing at one loop higher. These can be accounted systematically by replacing 
Eq. (4) with an effective action 5*;. 

Si = j[S{Vl(P) + Sm{Vl4>) - Sm{Vl4>)] . (5) 

Starting from this effective action the self energy S of the complete propagator can be 
calculated to any order in / and set / = 1 at the end of the calculation. The gap equation 
is obtained by demanding that S does not shift the mass m, i.e., 

\p2=.m2= (6) 

Furthermore, in order to get a real solution of m from Eq. (6) one has to ensure that the 
imaginary part of at = m? is zero. 

We apply the above mentioned ideas to (7^0^ theory (described by Eq. (1)) at finite 
temperature using real time formalism. Let us recall that in the real time formalism 
P, ^ the thermal propagator has a 2 x 2 matrix structure, the 1 — 1 component of which 
refers to the physical field, the 2 — 2 component to the corresponding ghost field, with 
the off-diagonal 1—2 and 2 — 1 components mixing them. The propagator used here is 
given by 



Dn{K) D,2{K) \^{ AoiK) + AfsiK) Ap{K) 
D2i{K) D22{K) ) [ Ap{K) Al{K) + Ap{K) 

where Aq{K) is the usual Feynman propagator at zero temperature 



(7) 



K'^ — + le 

Here Ap and are finite temperature corrections to the zero temperature propagator 
where 

Ap{K) = 2ix5{K'-m')nB{\m, A^(ir) = 27r5(ir^ - m^P ^B(|i^o|), (9) 

with the Bose-Einstein factor nB{Ko) = l/(e^''^°' — 1) (/? = 1/T). The complete self 
energy is given by the expression pi][ : 

ReS(p) = ReSn(p), ImS(p) = -—^-j^j:,,{p), (10) 

where En and S12 are the self energy of the 1 — 1 component and 1 — 2 component, 
respectively. 
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The graphs that are contributing to E(j9) up to two-loop order are depicted in Fig. 1. 
The contributions to ReS(p) coming from these graphs, although in a different context, 
have already been computed in Refs. |^, ^, |1T| . We find that the real part of self energy 
in the limit po = m and p ^ is given by 

f oo oo ^ 

ReE(po = m, p ^ 0) = -y^T' + g^T'l -ny + y^Y. Mv) - E Mv) 



n=l 



n=l 



+g T {2tt + aiy - 27iy \n y - a2y + a^y In y 
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(In y^f - (a4 + In y^ 



y n=l 



- (asy' - 2ny) ^ Bniy) + E E A^(?/)5„(y) 

n=l 771=1 n=l 

+ A[X{y) + I{y)]]+0{g') , (11) 



where y = f,g = ^,a, = 15.64-21.99 In (J), = 26.07-3.46 In (5)-2(ln ( J))2_i5,+ 
i52, as = 1.62 + i In (f ), = 2.54, a, = 0.98 + ln (J), AM = T;.[8'^'n\l + j^)-^ -y% 

5n(l/) = ^[(l + i.)-^-l], 

^1 = /o dxi dx2il - X2) In [(1 -X2 + jrz^^) - ^2(1 - ^2)], 

52 = /o' dx, rfX2(^ - 1) In [(1 -X2+ (Y^) - X2(l - X2)], 

^(y) = y' /r In {t + V^^), 

T(„.) _ „,2 TOO dt2 pp TOO dh 1„ fV^zlzVEEl] 

ny)-y Jl (e.t2_l)i^i^ Jl (eVH-l) in [^^^^^J. 

The only graph that will contribute to the ImS(p) is Fig. 1(h) and its contribution to 
Si2(p) is given as 

- iSi2(p) -- 



g^ f d^Ki d^K2 



6 J (27r)4 (27r) 



Ap{Kr)Ap{K2)Ap{p -Kr- K2] 



f d^ki d^k2 . UBiEkJuBiEk^] 



X U 



(po - Ek, - Ek,Y - E^ 



'p—ki — k2 



+ S (po - Ek, + Ek^f - El_,^_,^ 



+ 6 (po + Ek, - Ek^ f - El 



ki — k2 



+ s 



{po + Eki + Ek^y - El^f.^_f^,^ 



(12) 
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where = -^/k^^'+m^, = -^/k^^'+m^, and Ep_ki-k2 = y'(p - ki - ka)^ + m^. In 
the po — m and p ^ hmit, the expression for — iEi2(p), becomes, 

- iEi2(po = "^,p) = -——^ d^kid^k2-=r-E^nB{Ek^)nB{Ek^) 

(13) 

If ki makes an angle 6 with the direction then, the above expression becomes 

X ^5{x - xi)e{l - \xr\) + 6{x - X2)e{l - \x2\) 

+ 5{x - x^)9{l - Ixsl) + 5{x - x^)9{l - |a;4|)| (14) 

where 



X — cosO 



(m - Ek, - Ek,f - ki^ - W -m" (m - E^, + E^.f - W - W - 



^1 = 



2 



(m + Efc, - Ek,f - ki' - ka' - _ ^ (m + + EfcJ' - ki' - ka' - 
2fciA;2 '^"^ 2fciA;2 

Now we have to evaluate the range of integrations over ki and ka using this theta 
functions. 

1) 9{1 — \xi\) = 9{1 — xi)9{l + xi). Therefore xi < 1 implies 
^JEk^ -m^Ek^ - m[^J {E^^ + m){Ek^ + m) - ^J{Ek^ - m){Ek^ - m)) > 0, 
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and it is trivially satisfied for positive values of ki and k2- Similarly xi > —1 implies 

{Ek-^ - m){Ek2 -m) + kik2 > 0, 

and is also satisfied by both of ki and k2- Therefore, this theta function does not put any 
extra constraint on the integration range of ki and ^2- 

2) ^(1 - \x2\) = ^(1 - X2)e{l + X2). X2<1 implies 

{Ek, - m){Ek^ -m) + kik2 > 0, 
and it is trivially satisfied for positive values of both of ki and k2. Next X2 > —I implies 

k2{k2 - ki) > 0. 

Therefore it gives 9{1 — \x2\) = 6{k2 — ki). 

3) ^^(1 — l^al) = ^(1 — 2:3)^(1 + 2:3). Proceeding in the similar fashion one gets 

e{i-\x:i\) = e{k^-k2). 

4) e{l - \x4,\) = 9{1 - X4)9{1 + X4). X4<1 implies 

k2{k2 + ki) < 0, 

and this relation is not been satisfied for positive values of ki and k2. Therefore this theta 
function will not give any contribution to — iEi2(po = m, p — > 0). 
Using this theta function constraints we ultimately get, 

g'^ f°° dkik2 dk2k2 



•V f 9^ f^dkik2 f 



X 



24(27r)3io Ek, Jo E^, 

nB{E,,)nB{E,;)eP'''^J'eP^'^^'' 

|ns(Efc, + Ek^ - ^)e/3(^*i+^'=2-"»)/2 

+ e{k2 - k^)nB{Ek, - Ek, + m)e^(E,^-^'^^+my2 

+ e{ki - k2)nB{Ek, - Ek, + m)e^(^'=i-^'=2+'")/2| ^^g) 

Now if we interchange ki and k2 in the last term of the above expression it becomes 
identical to the second term, and we get. 



-zT,i2{po = m,p^0) -- 



24(27r)3 
+ 



[ Jo Ek^ Jo Ek^ 



dkoko 



/'^ ^e^^^.e^-/^n^(£;., - Ek, + m) \ 

ryk2 JO J^ki J 

xnB{Ek,)nB{Ek,) (17) 
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Finally the ImE(po = m, p — > 0) takes the following form: 



i Ql3m _ I 

ImE(po = m, p ^ 0) = - ^^^^^ ^uiPo = m, p ^ 0) 



—g'T'Hy) (18) 



where nsix) — l/{e^ — 1) and 



2 /* oo /* oo 

J{y) = 2 (ln{l-e-y)) -3 J dx\n{l - e~^) - 2 J dx{x - y)nB{x + y) (19) 

Therefore it is evident from Eqs. (11) and (18) that S(p) \(^ pq = m,p 0) is complex 
Consequently, the thermal mass(m) that we may obtain by solving the gap equation 
^(p) |( Po = ""^i P ^ 0) = has an imaginary part. However, it is clear from Eq. (10) 
that due to a multiplicative factor (e^^° — 1), ImS(p) is zero at Po = |p| and |p| — >■ limit 
and hence E(p) is real in this limit. 

At very high temperature y ^ 1 and the real and imaginary part of the self energy in 
this high temperature limit takes the following form: 



ReE(y) = -TV + 9'T'{^-7:y} + g'T'{3.7n'-'^--hy-hy' 

+ (— + 27^)W - hy^W] + 0(^^) (20) 

and 

ImE(y) = g'T'{^ + 7.9y - ^y' + "^yW - \{Wf} + 0(^^) (21) 

where, hi = 37.12 — oi, 62 = ^2 — 7.36 and 63 = 2.7 — 03. Therefore the gap equation 
in this limit is 



ReE(y) +iImE(7/) = 0. (22) 

Since y is complex, we set y^ — yR + iyj and substitute it in the high temperature 
approximated gap equation(22) to obtain the following two equations: 

yl-fR{yR,yi) (23) 
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yi = fi{yR.yi) (24) 

where ur and yi are the real and the imaginary part of the thermal mass respectively 
in the high temperature limit and the explicit form of and // are given in the appendix. 
The eqn.(23) and eqn.(24) are transcendental equations in two variables Ur and yi and 



we can solve this equation by the method of iteration |T^ assuming that our solution is 
correct up to order g^. We have started first with an approximate values of a pair of roots. 

After third iteration we find that the improved pair of roots are equal to the pair of 
roots obtained after second iteration. Therefore at very high temperature the real and 
the imaginary part of the thermal mass up to order is 



and 



m 



R 



+ 477 TTln T 
3 ' ^g 



(25) 



m, 



r27r^ vr /, /27r^\\2-| 47r^^, /1\ 47r^„/ /1\\2 , , 

[— - 3 )) 1^ ^ + ^"(^) - + ^'^^ 

We see that up to g'^ order both m|j and m\ are independent of the ultraviolet scale 
/i used in the theory. However from the structure of and //, it is evident that beyond 
g^ order they ought to be //-dependent. The result of m|. up to order g^ matches with 
that obtained by Parwani[|l^], however his g^ order term is /x dependent. As a result the 



two loop real thermal mass obtained in |^ becomes unstable below some characteristics 
scale of the order of g^^'^T . 

In this letter we have studied the finite temperature gap equation in massless g^(^'^ 
theory and its nature of solutions. We find that the thermal mass up to two-loop order 
which one may obtain self-consistently by solving this equation is complex. In the massless 
(7^0"^ theory the typical one-loop thermal mass is of order gT and the expected contribution 
to the mass from the two loop level would be of order g^T . It is worth mentioning that 
there may be some non-pertubative features that really start from this g^T scale which 
makes the computation unreliable beyond g^ order. It will be also interesting to extend 
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this method to 3 + 1 dimensional QCD where the generation of magnetic mass is quite 
problematic due to IR divergences in the two-loop level. 

I gratefully acknowledge Prof. S. Mallik for helpful discussions. I also thank the referee 
for pointing out a major mistake in the evaluation of the imaginary part of the self energy 
and also for his suggestion to solve the gap equation. 
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Appendix 



fniyR: yi) 



+ 



yl 



hi 



1 



^\/yR + yi + ^yl 



^2?/^ + (tt + ^)ln(z/^ + yj) - ^hylhi(y\ + y^) + 63Z//tan-i( 



-7.9 



rl 



1 



TTrl 



^^yi + y^- ^yl 
\\/yR + yi + lyR 



^yl 



-tan — 

3 ^yl 



^^y'R + y^)} 



(A.l) 



fi{yR,yi) = -7r[^v^^Tyf - '} 



+4^ + 



27r3 


p 1 


yl 1 


3 


^^^/yi + y^ 





-hx\-\jyR^y\- -yV\ ' - 



+ (27r + ^)tan-^(4) - ^?^3y|ln(|/A + y^) - h,yli^^-\ 



^yR 



yl 



TT 



(A.2) 
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Figure caption 

Fig. 1. Graphs contributing to the two point function. Sohd cross and sohd dot 
denote vertices for renormaUsation counterterms of 0^ and 0^ type respectively, while a 
dashed cross represents the additional 0^ insertion. 
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(b) 




(h) 
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